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Abstract. The yields of heavy flavored hadrons emitted from strangeness- rich 
QGP arc evaluated within chemical non-equilibrium statistical hadronization model, 
conserving strangeness, charm, and entropy yields at hadronization. 



1. Introduction 

A relatively large number of hadrons containing charm and bottom quarks are expected 
to be produced in heavy ion (AA) collisions at the Large Hadrons Collider (LHC). We are 
interested in how the high strangeness yield influences heavy and multi-heavy hadron 
yields (containing more than one heavy quark) [I]. Our work is differing from other 
recent studies [21 [3] , since these assume that the hadron yields after hadronization are 
in chemical equilibrium. However, we form the charm c (and bottom b) hadron yields 
in the statistical hadronization approach based on a given abundance of u, d, s quarks 
fixed by the bulk properties of a chemically equilibrated QGP phase. Since at LHC, 
and also RHIC, conservation of baryon number is not a vital element of the analysis, 
only entropy is preserved in the following, considering u, d, u, d yields. 

Our approach is justified by the expectation that, in a fast break-up of the QGP 
formed at RHIC and LHC, the phase entropy and strangeness will be nearly conserved 
during the process of hadronization. If and when the volume doesn't change appreciably 
during hadronization, the entropy conservation determines values of hadron yields after 
hadronization. For QGP in chemical equilibrium, this implies that the final hadron 
yields are in general not in a chemical equilibrium. In this work, we will show how 
this influences the relative yields of heavy flavored particles in the final state. We 
will investigate in quantitative terms how such chemical non-equilibrium yields, in the 
conditions we explore well above the chemical equilibrium abundance, influence the 
expected yields of single, and multi-heavy flavor hadrons. We show that there are 
profound differences in hadron yields between the assumed hadron equilibrium and our 
approach. 
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2. Statistical hadronization model with conserved yields 

During a fast transition between QGP (Q) and HG (H) phases strange i — s and 
heavier i — c and i = b quark flavor yields are preserved, as is the specific, per rapidity, 
hadronization volume, and in scaling limit, we also preserve the entropy per unit of 
rapidity 

dNf__dNl_ dS^__dS^^ dV^__dV^_ 
dy dy ' dy dy dy dy 

From now on we omit the upper index H. All parameters are for hadron side, unless 
there is an upper index Q. 

The number of particles of type l V with mass per unit of rapidity is, in our 
approach, given by: 

dNi dV 

~dy~ = Ti ^ n * • (2) 

Here, T, = X^i is the particle i yield fugacity, dV/dy is the system volume associated 
with the unit of rapidity, and nf 1 is a Boltzmann particles phase space density, 

nT = 9i J exp(- v /p 2 + m 2 jT) = g i ^^K 2 (m i /T), (3) 

where gi is the degeneracy factor, T is the temperature and p is the momentum. 

The particle fugacity A, = e Mi//T is associated with a conserved quantum number, of 
strangeness, charm, bottom, and baryon number: i = s, c, b. Antiparticles have negative 
chemical potential and thus A = A -1 . The A evolution during the reaction process is 
related to the changes in densities due to dynamics such as expansion. The phase space 
occupancy 7, is the same for particles and antiparticles. Its value can change as a 
function of time even if the system does not expand, for it describes buildup of the 
particular particle species. For this reason ji very often changes rapidly, while A is more 
constant. It is 7^ which is most sensitive to the time history of the reaction, except if 
full chemical equilibrium can be established, in which case 7$ — ► 1. 

The yields of hadrons after hadronization are given by Eq. (|2j) , the unknown phase 
occupancies on hadron side for strangeness 7^ and heavy quarks j c ^) can be determined 
comparing quark yields in Q, H phases. The 7 s /7 9 ratio depends mostly on strangeness 
to entropy ratio s/S after hadronization, see figiU obtained using SHAREv.2.0 [5]. For 
LHC the expected ratio s/S = 0.038 [6j at T = 140-180 MeV which implies in the 
hadron phase 7^/75 = 1.8-2 [7]. 



3. Yields of heavy flavored hadrons 

3.1. D, Ds, B, Bs meson yields 

The predicted D and D s mesons yield ratio D/D s is shown in figure [2] on left, as a 
function of ratio ^s/lq for T = 140,160,180 MeV. A deviation of 7^/7^ from unity 
(from chemical equilibrium) leads to a noticeable change in the ratio D/D s , which is 
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Figure 1. (color on line) Strangeness to entropy ratio, s/S, as a function of 7s / 7 g- 
(solid line, blue) for T = 200 MeV, S H = S Q -> j q = 0.083; (dashed line, blue) 
for T = 170 MeV, S H = S Q -» j q = 1.15; (dash-dotted line, green) for T = 140 
MeV, S H = S Q -» 7 g = 1.6; (dot marked solid, violet) for 7g = 1,T = 170. 

proportional to the inverse of 7 s /7 g ratio. Ratio 7 s /7 g grows with growth of strangeness 
to entropy ratio s/S. For LHC, s/S may reach a high value, and if s/S = 0.04 then 
the ratio 7 s /7 g is increased by about a factor 2 and D/D s ratio is decreased by a factor 
2 compared to the chemical equilibrium model. For B, B s mesons the results are the 
same as for D, D s mesons. 

3.2. Heavy baryon yields 

On the right, in figure EJ we show ratios cqq/cqs = (A c + S c )/S c and cqq/css = (A c + 
S c )/fi c as a function of 7 s /7 g for T = 200 MeV (dash-dot line), T = 170 MeV (solid line) 
and T = 140 MeV (dashed line). The ratio cqq/cqs is inversely proportional to 7 s /7 9 
ratio, similar as for D s /D ratio. However, cqq/css oc (7 s /7 g )~ 2 , further enhancing the 
effect of chemical non-equilibrium in QGP hadronization, with multistrange-charmed 
baryons produced at similar yield as nonstrange-charmed baryons. 

3.3. Yields of hadrons with two heavy quarks 

The yields of hadrons with several heavy quarks somewhat depend on hadronization 
condition, in particular on the reaction volume and on the total assumed charm yields. 
If we normalized yield of, for example, hidden charm mesons on their yield only 
slightly depends on total charm multiplicity. The yields of hadrons with two heavy 
quarks are approximately proportional to \/{dV/dy) [Sj, because 7^ for heavy quarks is 
proportional to l/(dV/dy) [8], and only one power is cancelled by the proportionality to 
volume. The results we present are based on QGP reference state with dV/dy = 800 fm 3 
at T = 200 MeV, for s/S = 0.04 and dV/dy = 600 fm 3 at T = 200 MeV, for 
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Figure 2. (color on line) Left: D/D s ratio as a function of 7^ /7g f° r T = 140 
MeV (blue, dashed line), T = 160 MeV (green, solid line) and T = 180 MeV 
(red, dash-dot line). Right: The ratios cqq/cqs = (A c + S c )/H c (lower lines) 
and cqq/css = (A c + S c )/Q c (upper lines) for T = 200 MeV (dash-dot line), 
T = 170 MeV (solid line) and T = 140 MeV (dashed line). 



s/S = 0.03. Volume expands with drop in temperature preserving dS/dy. This 
corresponds considering entropy to a total final hadron particles multiplicity of about 
5000 for s/S = 0.04. 

In figure [3] we show the yield of hidden charm mesons cc (sum over all states of cc) 
normalized by the square of charm multiplicity per rapidity N% (N c = 10 for LHC and 
N c = 3 for RHIC) as a function of hadronization temperature T. We consider the cases 
s/S — 0.03 and s/S — 0.04 for LHC. We see that while the chemical equilibrium model 
predicts a yield falling with hadronization temperature T, the fixed s/S yield remains 
roughly independent of T. Remarkably, we see comparing right and left sides of figure [3] 
that in our model the expected yield of cc mesons decreases with increasing specific 
strangeness and light quarks yields. The multiplicity of light quarks is also larger than 
in chemical equilibrium for T < 200 MeV in the order to secure entropy conservation 
during hadronization. This new mechanism of charmonium suppression occurs due to 
competition with the yields of D and D s . D and D s in effect deplete the pool of available 
charmed quark pairs, and fewer hidden charm cc mesons are formed. For larger T, right 
part of figure El the multiplicity of light quarks becomes smaller than equilibrium. This 
results in an enhancement of cc yield in comparison with chemical equilibrium. 

In figure H] we quantify a more systematically this result. We compare the J/9 yield 
to the chemical equilibrium yield (J/9)/ (J/9 eq ), as a function of js/jq, each line is at a 
fixed value 7 ? . This ratio is: 

J/* = N hld = jl_ 
J/9 eq N hideq 7c 2 ' 
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Figure 3. cc/N^ relative yields as a function of hadronization temperature 
T. Dashed lines are for chemical equilibrium, solid lines for prescribed s/S in 
QGP. Left panel: dV/dy = 600 fm 3 at T = 200 MeV, with s/S = 0.03 for the 
nonequilibrium curve. Right panel: dV/dy = 800 fm 3 at T = 200 MeV with 
s/S = 0.04. 



(Jjty)/(Jjty cq ) always decreases when ■y s /lq increases, as we expect. We see that when 
we have small yields of light q, s quarks, we of course find J/$>/ J/$ eq > 1. This happens 
e.g. when hadronization is at T = 200 MeV. In all other and 'reasonable' hadronization 
scenarios we find suppression (J/4 r )/(J/4 r eq ) < 1. 



4. Conclusions 



We have considered the abundances of heavy flavor hadrons within the chemical non- 
equilibrium statistical hadronization model. We studied how the (relative) yields of 
strange and non-strange charmed mesons vary with a QGP produced fixed strangeness 
content. A considerable shift of the yield from non-strange D(B) to the strange D S (B S ) 
and similar for strange-charmed baryons is expected for LHC since s/S may reach value 
about 0.04, which upon hadronization means ^s/lq — 2. 

Another important consequence of this result is that we find a relative suppression 
of the multi-heavy hadrons, including J/^f, for hadronization temperature lower than 
T = 200 MeV, because strangeness and light quarks are above chemical equilibrium 
in HG after hadronization. We compared the yields to the expectations based on 
chemical equilibrium yields of light and strange quark pairs and found that the J/^> 
is suppressed, the mechanism being that abundantly produced light quark-charmed and 
strange- charmed hadrons deplete the pool of available charmed quarks. 
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Figure 4. (Color on line) Ratio J/3/ / Jfl> eq = Jc/^ceq as a function of ^ s /lq at 
fixed value of j q and if required, entropy conservation. Shown are: T = 200 
MeV at 7 g = 0.83 (dot-dash line, blue); T = 170 MeV at -y q = 1 (dashed line, 
green) and at j q = 1.15, (solid line, red); at T = 140 MeV, j q = 1.6 (solid 
dotted line, purple) 
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